Entanglement properties of quantum spin chains 
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We investigate the entanglement properties of a finite size 1+1 dimensional Ising spin chain, and 
show how these properties scale and can be utilized to reconstruct the ground state wave function. 
Even at the critical point, few terms in a Schmidt decomposition contribute to the exact ground state, 
and to physical properties such as the entropy. Nevertheless the entanglement here is prominent 
due to the lower-lying states in the Schmidt decomposition. 



PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Pq, 05.70. Jk 



O 
O 

(N 

^ ■ 
< 

o 

cn 

> 
m 
m 
(N 
(N 
O 

O 

:^ 



I. INTRODUCTION 

The EPR argument ,1] and the Bell inequalities 
were, albeit with almost 30 years in between, acknowl- 
edgments that quantum theory exhibits the strange cor- 
relations known as entanglement. However, only in 
the latest few years has it been known that entangle- 
ment is a resource that can be utilized in quantum 
computing, 3'] and is thus a central subject in the con- 
tinually expanding field of quantum information theory. 
Furthermore, entanglement has been shown to be a fun- 
damental feature in quantum phase transitions, some- 
thing that has spawned a whole new field of research 
IllSliaSlllillSlilliilS- This article will focus 
on entanglement in both quantum critical and non criti- 
cal systems, the structure of entanglement in condensed 
matter systems being at its most complex in critical sys- 
tems. 

The determination of the ground state in a quantum 
system such as ID spin chains is highly complex, but the 
last years have seen the development of techniques such 
as DM RG IT^ fl^ and the recent entanglement renormal- 
ization |15| which solve the problem efficiently. In both 
cases the entanglement in the problem seems to be the 
key to the methods' success over traditional renormaliza- 
tion schemes. Moreover, for a pure state, any long range 
correlations require the existence of entanglement in the 
state. 

A quantum system with wave function j^*) in a Hilbert 
space H, which is partitioned into two subspaces 
Ti.sg, can be written as a Schmidt decomposition 



I*) 



(1) 



where \ip^^) € Ti^ and vice versa. The coefficients A„ 
are real, positive c-numbers and the states are mutually 
orthogonal, (V'^IV'm) = ^nm- The upper limit x> the 
Schmidt number, is a brute measure of the entanglement 
between the two subsystems, running from one (no entan- 
glement), maximally reaching the dimensionality of the 



smallest of the two Hilbert spaces. The reduced density 
matrix of (say) system can then be written 



= Tr^|vI')(M/| = ^A„|V'r))(V'i' 



and equivalently for £§. Hence, knowing the Schmidt de- 
composition of the wave function is equivalent to knowing 
the basis in which the reduced density matrix is diago- 
nal, modulo phases. We will in this article show that the 
effective Schmidt number Xoff is much smaller than the 
possible maximum, and thus the number of terms con- 
tributing significantly to the wave function is surprisingly 
small. Remarkably, this also is true at highly entangled 
points in the phase space, such as at criticality. 

The entropy of the wave function |^) can be measured 
as the entropy of the subsystem. 



5^ = -Trp(^)log2p(^)=-5]A. 



log2 A„, 



(2) 



increasing from zero if the two subspaces are entangled. 
Hence, computing the eigenvalues A„ of the reduced den- 
sity matrix is vital to the understanding of entangled 
states. 



II. ISING MODEL 

The quantum Ising chain in external field h in 1+1 di- 
mension is a good benchmark tool for the analysis of en- 
tanglement since its properties are extremely well known. 
The model is defined by the Hamiltonian 



Ising — 



N 

E 

ri=l 



(3) 
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on N lattice sites. We operate with open boundary 
conditions (OBC), since this will ease the formalities 
when partitioning the system. This will weaken the 
phase transition in the sense that the conformal sym- 
metry is broken compared to periodic boundary con- 
ditions (PBC) |l6j, but the phase transition in the 
thermodynamic limit will prevail. Also, Calabrese and 



2 



Cardy have found that one can compute conformal sig- 
natures in the OBC case provided the boundary condi- 
tions are conformal The model has a phase tran- 
sition at /i = 1 between the product state \'^oo) — 
I It • • • t ) for h > 1 and the Schrodinger cat state 
1*0) = ^ (I • • • ^) + I • • • ^)) when h<l. 
Here |t) is the eigenstate of cr^ with eigenvalue 1, and 
1^) are eigenstates of with eigenvalues ±1. Note that 
the transition is between unit entropy in the low field 
limit and zero in the high field limit. 

Doing a Jordan- Wigner transform as sketched in Refs. 

llSj , we can map the model onto a string of non in- 
teracting fermions, and thus compute the eigenvalues for 
the Ising model (and a larger class of models) in what 
resembles the thermodynamic limit, typically a few hun- 
dred particles. In essence, we define the N fermionic 
operators 




(4) 



and the 2N Majorana fermions 
1 ... 



1 



These fulfill {7m, In} = Smm and are delocalized in terms 
of the original lattice of fermions. The Majorana fermions 
diagonalize the Hamiltonian in the sense that 



H 



with an Hermitian 2N x 2N matrix C. Next, define 
the imaginary and anti symmetric correlation matrix 
Tij = ([7i,7j])- We will consider only the ground state, 
and thus the expectation values are to be taken in the 
ground state. In this sense, our model is defined by the 
matrix C while the state is defined in the matrix F. The 
Majorana fermions have a two-to-one correspondence to 
the fermions iQl, and hence tracing out a particle from 
the system amounts to removing the two adjacent rows 
and columns in F corresponding to the particle. Trac- 
ing out N ^ N' particles this way, we end up with the 
2N' X 2N' matrix Fy. This latter matrix can be made 
block diagonal through an orthogonal transformation O, 
such that we can define new Majorana fermions 7„ that 
fulfill 



N' 



OTfO = {([7„7,])} = 



fc=i 











(5) 



with < S,k < 1. The transformation O is the same trans- 
formation that block diagonalizes the Hamiltonian ma- 
trix C. The matrix F corresponds to the state described 
by the reduced density matrix p' , with the N — N' parti- 
cles traced out. The fermion operators o„ corresponding 
to 7„ will diagonalize the reduced density matrix such 



that the eigenvalues thereof can be computed. The eigen- 
values are determined by the set of binary occupation 
numbers rj = {uk}, fc = 1, 2, . . . , N' , and nk = {0, 1}, \ri) 
being an eigenstate of afca|, with eigenvalue Uk- Thus \rj) 
is also an eigenstate of p' and the reduced density matrix 
is p' = J2ri W''l){v\ with eigenvalues 



A. 



N' 

n 

fe=i 



(6) 



Finally, the entropy becomes 
^ /I 



k=l 



(7) 



where H(x) = —x logj x—{l — x) log2 (1 — x) is the binary 
entropy function. 

The entropy of the ground state as measured by Eq. 
(0) when the system is partitioned into two equal parts, is 
shown in Fig. ^ The increased entropy around the criti- 
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FIG. 1: The entropy in the quantum Ising model as function 
of the external magnetic field h for = 10 (lower line) and 
= 100 particles. In both cases half system is traced out, 
from the edge. The letters refer to the places of investigation 
in Fig. 121 The critical point (in the thermodynamic limit) is 
shown as a vertical line. 

cal point is a hallmark of the quantum phase transition, 
though this does not show in the small system. 

Figure |21 shows the magnitude of the decreasingly or- 
dered eigenvalues A^. Note that there are 2^/"^ = 128 
possible eigenvalue contributions, while this shows that 
only a very few of these contribute to the wave function, 
and thus to the entanglement as measured by the entropy 
S. 

The main eigenvalues at critical point decrease roughly 
exponentially, as A^ ~ e~^", while off criticality only the 
very first contribute significantly. In the low field limit, 
the highest eigenvalues are paired due to the Schrodinger 
cat nature of the state, while in the high field limit there 
is one single main contribution, which is the product state 
l^'cx;)- With PBC, the eigenvalues would be paired also 
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FIG. 2: The decreasingly ordered eigenvalues of the re- 
duced density matrix of the Ising model's ground state when 
tracing out half size. Data for A'^ = 50. The upper line (•) 
denotes the critical values when h = he — 1, a.t (B) in fig- 
ure (A) denotes the Scrodinger cat state with h = 0.5, 
at (A); and (a) denotes the approximate product state with 
ft = 2, at (C). The fitting line to the critical values ~ e~^" 
is shown. Only the 10 largest eigenvalues are shown since 
numerical errors become prominent after this. 



for h > 1, except for the single main eigenvalue due to 
the translational symmetry. This pairing also occurs in a 
conformal bosonic chain 1(11, but is not a general prop- 
erty of a non critical system. 

Hence we can approximate the wave function by re- 
stricting the sum Q to some upper limit x' < % to make 
the new wave function 
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where c is the central charge of the conformal field theory 
corresponding to the phase transition. In the Ising case 
c = 1/2. 
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FIG. 3: From the top, the figures show the four smallest 
eigenvalues of the correlation matrix in Eq. the four 

largest eigenvalues of the density matrix as computed from 
Eq. @, and the corresponding entropy contributions s„ = 
— A„ log A„ . For the two lower graphs, equal point types refer 
to the same eigenvalues. System is at the critical Ising point 
and is traced at half size. 



where the new coefficients are determined by normal- 
ization, = A„/(l — e). We define an error as 



e = E A„ 

n=x' + l 



which measures the difference in the eigenvalue sum of 
the two wave functions. The overlap becomes 



Figure 13 shows how the largest eigenvalues scale at the 
critical Ising point. The largest eigenvalue A{(A^) de- 
creases with system size at criticality, while the other 
increases. The largest decrease roughly as A];(A^) ~ 
— 0.027 logiV, while the second increases roughly as 
A^(iV) - 0.021 log TV, both within the range shown in 
Fig. |31 However, all entropy contributions increase with 
N, thus contributing to the overall divergence of the en- 
tanglement entropy. 



III. SCALING OF EIGENVALUES 

We investigate how the eigenvalues scale with increas- 
ing system size. Off the critical point there is no scal- 
ing with system size of the eigenvalues since the entropy 
saturates at some value. However, on, or near the criti- 
cal point the entanglement entropy diverges as predicted 
by conformal field theory. Given that we trace out a 
constant fraction of the entire system, the entanglement 
entropy diverges up to an additive constant as [a. 11^ 



IV. APPROXIMATING THE GROUND STATE 

There exists successful techniques to approximate the 
ground state of a quantum s yste m, the most prominent of 
which is the DMRG scheme [l^ . Rather than taking only 
ground states of a sub lattice when renormalizing the sys- 
tem, this method takes into account possible entangled 
states as well, producing unprecedented accuracy. The 
role of entanglement in DMRG is still under investiga- 
tion, but it seems clear that it is vital to the success of 
the scheme [T^ . The results of our work here are in- 
dicative that the main contributions to the ground state 
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are indeed entangled at the critical point, along with few 
of these terms contributing to the actual wave function. 
These results are not directly applicable to an improve- 
ment of the DMRG algorithm, but rather indications as 
to the success of DMRG. Furthermore, any prospective 
technique to find the ground state of a quantum system, 
even in those areas where known techniques fail, needs to 
understand the nature of the wave function in the system, 
and we believe that the entanglement properties would 
play an important part in such a method. In particu- 
lar, the entanglement properties is precisely what distin- 
guishes a quantum many-body system from the classical 
counterpart, and entanglement must therefore be an es- 
sential part of any such method. 

Having found that only few terms contribute to the 
entanglement in the ground state, we query how well we 
could possibly approximate the ground state of the full 
system by the first few terms in the Schmidt decomposi- 
tion. To this end, we focus on the few-particle case, where 
the exact wave functions can be computed explicitly. We 
split the open spin chain in two partitions, and compute 
the overlap between the actual wave function and the 
Schmidt expansion of the two subsystems, ordered de- 
creasingly on those terms with the largest overlaps. 

A priori it is clear that in the zero field limit Oi = O2 = 
l/\/2, while any higher terms vanish. In the high-field 
limit, Oi — 1 as the only non vanishing overlap. We 
assume that the wave functions \ipn^) are ordered with 
decreasing Schmidt number, and the results are shown 
in Fig. H 




FIG. 4: The four largest overlaps On, n = 1, 2, 3, 4 (from top 
to bottom) with magnetic field h and system sizes A'' = 10 
(full lines) and A*' — 100 (dashed lines). Again, the critical 
point is indicated with the vertial line. 

Since the fermions Ofc corresponding to the Majorana 
fermions defined in Eq. jSJ diagonalize the reduced den- 
sity matrix, the overlap follows straightforwardly. The 
probability of the state jyy) given the density matrix 
p' = A^|7y)(?7| is A^, and the overlap between this and 
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FIG. 5: The three main entropy contributions s„, n = 1,2,3 
(from top to bottom) in the Ising model with magnetic field h. 
Full lines are for A' = 10, while dashed lines are for A' — 100. 



the ground state will therefore be ^/A^, or 

k 



Most prominently, the Oi overlap increases monotoni- 
cally as h increases. That is, even at criticality this first 
term consisting of the ground states at each side approx- 
imates the actual wave function even better that it does 
for the cat-state where the overlap naturally is 1/V2- 
The increased entropy at critical points is mainly due to 
other terms rising around criticality, terms that are not 
part of the entire wave function at non critical points in 
the phase space. 

Investigating the Schmidt coefficients' entropy contri- 
butions, the results are shown in Fig. [S] Hence we see 
that apart from the two major contributions that are 
equally prominent bX h ^ Q, there is one more contribut- 
ing at criticality. The effect is mainly pronounced with 
big systems, but a trace is also seen at = 10. The 
increased entropy at the critical point is contributed to 
by the largest term, but also by the appearance of terms 
that are zero at all noncritical points. 

The error in the overlap and in the entropy that arises 
from only choosing x' < X terms in the Schmidt decom- 
position is zero (or very close) off criticality, since the 
Schmidt rank is very small. However, on criticality x in 
principle extends to the dimension of the Hilbert space. 
However, the errors in the overlap 50 = 1 — X]n=i 

and entropy SS ^ S — J2n=i plotted in FigEl We 

see that the errors grow linearly with N when N > 100. 
The number of Schmidt terms needed to compute the 
entropy up to an error of lO"**, or the effective Schmidt 
number is only 4 up to a system of iV = 400 sites. 
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FIG. 6: The errors SO (o) and 5S (•) as defined in the text 
with different system sizes at the eritieal Ising point. For 50 
we have chosen x' — 4 and for 55 x' = 3. In all cases the 
system is traced at half-size. The lines for A'^ > 100 are linears 
JO ~ 5 X IQ-^iV - 3.8 X 10"* and 5S ~ 1.75 x 10"® AT - 8.1 x 
10~^. For small systems the linearity does not hold. 

Conclusions 

We have seen that the complex structure of a wave 

function at critical point in the Ising model conies from 
the emergence of several terms in the Schmidt decompo- 
sition of the wave function. Nevertheless, remarkably few 
terms actually contribute to the expansion, and physical 
properties such as the entropy can be extracted using 



very few terms, even at criticality. Mainly, the rapid 
decrease of the Schmidt coefficients on or off criticality 
identifies this effect. We have detailed the effects with 
both large and small systems, and the primary conclu- 
sions hold in both cases. 

The matter that these few eigenstates contribute to 
physical properties such as the entropy and also the over- 
lap with the true ground state signifies the fact that find- 
ing the true ground state, or an approximation to it, is 
not a complicated task per se, we known that only a few 
of the eigenstates of the Hamiltonian would suffice for 
a sufficiently good description of the state. The results 
in this paper also indicate that these states are indeed 
entangled, and thus any algorithm that does not include 
entanglement is bound to fail. 

The effective Schmidt number that is needed to find 
a state within some limit is small, on the order of 10 
even up to several hundred particles. Thus the effec- 
tive Hilbert space is dramatically reduced. However, we 
have not addressed the problem of determining the un- 
known ground state of a system, only pointed out that 
the problem is not consistently hard given the correct (so 
far unknown) approach. 
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